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equations with time-dependent dissipation 



Jens Wirth*t 



Abstract 

We consider the initial- value problem of abstract wave equations with weak dissi- 
pation. 

We show that under conditions on the dissipation coefficient and its derivative the 
solutions to the abstract dissipative equation are closely related to solutions of the free 
problem multiplied by a decay function. This paper gives the counterpart to a recent 
paper of T.Yamazaki [Adv. Differential Equ., ll(4):419-456, 2006], where effective 
dissipation terms and the relation to the corresponding abstract parabolic problem are 
considered. 

AMS Subject Classification: 35L90, 35L15 

1 Introduction 

Let H be a separable Hilbert space with norm || • || . Let further A be a closed and self-adjoint 
non-negative operator with domain T>(A) C H. Note, that the domain of A becomes itself 
a Hilbert space if we endow it with the graph norm = 1 1 xi 1 1 2 + ||^4u|| 2 ~ ||(^4)u|| 2 . To 
simplify notation we set for the following A = \[~A. 

We consider for a positive valued C 1 -function b = b(t) on [0, oo) the abstract dissipative 
wave equation 

\ u" + 2b(t)v! + Au = 



(i- 1 ) 

|u(0) =ui& 23(A), it'(0) = u 2 e H 
and compare its solutions to the corresponding free problem 

|v(0) = «i e 23(A), u ; (0) = v 2 E H 

in the abstract energy space associated to A under decay assumptions on the coefficient 
function b. 
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The first result is as follows. If we assume that b G i 1 [0, oo), then the solutions of (jl.ip 
are asymptotically free in the energy norm. This means, if we set E to be the closure of 
(D(A) Ker A) x H with respect to the norm 

\\{ui,u 2 )\\ E = \\(Au 1 ,u 2 )\\hxH, (1-3) 
the following theorem is valid: 

Theorem 1.1. Assume b G L [0, oo). T/ien i/iere exists an invertible operator W+ € C(E) 
of the abstract energy space, such that for (m,u 2 ) G E and (vi,v 2 ) = W+(ux,U2) the 
corresponding solutions to (jl.ip and (|1.2j) satisfy 

/CO 
6(r)dr - (1.4) 

as t — > oo. 

Remark 1. Note that the solutions to the free problem have the conservation of energy 
property, i.e. it holds for all t G M that = ||(fi, U2)[|.E- Thus the statement 

of Theorem 11.11 implies the non-decay to zero of the energy of solutions to the perturbed 
problem (II. lh together with the existence of asymptotically equivalent solutions to both 
problems. 

Remark 2. Note, that the operator W + in this theorem is the inverse of the M0ller wave op- 
erator. In this point our notation differs from [3]. This inverse operator is more appropriate 
to describe asymptotic properties of solutions to problems with perturbed coefficients. 

The aim of this article is to generalise this result to the case of so-called non-effective 
weak dissipation terms b = b(t) with 

(Al) b(t) > for all t G [0,oo), 

(A2) \b{t)\ < Ci(t)- 1 and \b'(t)\ < C 2 {t)- 2 , 

including non-integrable coefficients b £ L 1 [0,oo). In this case the abstract energy \\(u, u')\\e 
of the solution decays to zero. We will show that if is not an eigenvalue of A then the 
solutions can be written asymptotically as product of a free solution and a time-dependent 
function describing the energy decay rate. To be more precise, let us formulate a theorem 
(which is contained in Theorem 13. ID . 

Theorem 1.2. Assume (Al), (A2) together with lim sup^^ tb(t) < i and Ker A = {0}. 
Then to u\ G T>(A) and u 2 G H there exist corresponding data (vi,v 2 ) = W + {u\,u 2 ) G E 
such that the corresponding solutions to (jl.ip and (jl.2p satisfy 

||A(t)(u,u')-(i;y)||i5-0 (1.5) 

as t — > oo, where the auxiliary function X(t) is given by X(t) = exp(J Q * 6(r)dr). 
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Remark 3. If b L 1 [0, cxd) the integral in equation (jl.4p becomes infinite. In order to 
obtain a relation between the solutions u and w we have to pay. This can be seen from 
the differences of the statements. The operator W+ maps T>(A) x H continuously to E, 
it is in general not bounded from E to E. Furthermore, we cannot give a uniform rate 
of asymptotic equivalence. The auxiliary function \{t) describes the decay of energy. If 
Ker A = {0}, we get the two-sided estimate 

ll(«>«0ll*~ j^- (1-6) 

Since the convergence rate in Theorem 11.21 is not uniform in the data, the constants in this 
estimate depend in a (non- linear) way on the data (u\,U2) and do not yield a two-sided 
norm estimate for a corresponding solution operator. 

Example 1. One typical class of coefficient functions b = b(t) for our approach is 

b(t) = ^ , /x>0, n = 0,l,... (1.7) 

(1 + t) log(e + *)••• log N (eN + t) 

with iterated logarithms log^ = log, logf rt+ ^ = logolog^ and corresponding iterated 
exponentials = 1, e^- n+1 ^ = e^'™ 1 ). In this case X(t) ~ ^log' n l(e[ n l + t^j . The decay rate 
for the abstract energy becomes arbitrary small in the scale of iterated logarithms. □ 

Recently, T.Yamazaki, [11], |12j . considered the related problem in the case of effective 
weak dissipation, i.e. if tb(t) — * oo as t —* oo together with b(t) — * and suitable conditions 
on the first derivative. In this case, a relation to the free problem (|1.2p does not occur. 
Instead , there arises a close relation to the abstract parabolic equation b(t)w' + Aw = 0. 
While the purpose of [TT], [12] was to demonstrate the close relation of abstract wave 
equations with effective dissipation and the corresponding abstract parabolic problem, our 
purpose is to show what happens if the influence of the dissipation becomes less strong and 
how the abstract parabolic type asymptotics changes to an abstract wave type asymptotics. 
The used methods are based on [2]. 

This paper is organised as follows. First we will sketch one proof of Theorem 11.11 and 
give some remarks on possible generalisations. Then in Section [2] we will summarise the 
construction of representations of solutions to (jl.ip by means of a spectral resolution of the 
operator A and discuss its consequences for the decay rates of the energy. In Section [3] we 
state the main result of this paper and in Section H] we discuss some possible applications 
of the abstract results in concrete settings. 

The main part of this paper follows the author's treatment from [9] combined with 
ideas from [2] . In that paper precise L p —L q decay estimates are proven for solutions to the 
Cauchy problem of a wave equation with non-effective time-dependent dissipation. The 
approach was based on a diagonalisation scheme applied to the full symbol of the operator. 
Here we formulate the Cauchy problem in an abstract setting and base our representation 
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of solutions on a spectral theorem for the operator A. Besides this, there are two essential 
differences to the treatment in [9] . On the one hand the derivation of dispersive estimates 
is based on stronger assumptions on the coefficient function (we need estimates for n + 2 
derivatives of the coefficient function for dispersive estimates) and requires more steps of 
diagonalisation, while in the Hilbert space setting estimates for b and its first derivative 
are sufficient. Second difference to [9] is that we obtain the modified scattering result not 
only for lim sup^^ tb(t) < 1/2. We use a new idea to impose conditions on the data in 
order to obtain a similar result for lim inf^oo tb(t) > 1/2. 

1.1 Outline of the proof of the classical scattering result 

Theorem 11.11 can be reduced to abstract results already known from the literature. Using 
a substitution of the time variable t = t(r) = a(s)ds such that logA(i) = Jq b(s)ds = 
log a(r) reduces (jl.ip to the abstract problem w" + a 2 (r)Aw = 0. Because we have cl 2 {t) — > 
J* °° b{t)dt < oo and f£°(a(T) — a(oo))dr = J °° r(t)b(t)dt < do, the technique of A.Arosio, 
[I], can be applied. Transformation back yields Theorem II .li 

For completeness and because the idea of the proof is basic for our treatment in sections 
[2] and [3] we give an independent proof of Theorem II .11 

Proof of Theorem \l.l\ We use the canonic identification J : E — (H Ker A) x H of 
the energy space with cli?(A) x H and write the Cauchy problems in system form. This 
yields 

A f Au \ _ f Au '\ _ ( A \ f Au ^ 

dt\u'J~ \u"J~ \-A -2b(t))\u' 

We denote by £ (t, s) and £o(t — s) the corresponding semi-groups of operators, 



(li 



3^b(*)=(_ A A )«foW, Sq(0) = I e C(H x H), (1.9) 

^£(t, s) = ( £(t, s), £(s, a) = l€ C{H x H). (1.10) 



dt K ' ' \-A -2b(t) 

The free semi-group £o(t) is unitary (because its generator is skew); while under assumption 
(Al) the semi-group £(t, s) is contractive^ If we consider the operator 

Q(t, s) = £ (s - t)£{t, s) G C(H x H), (1.11) 

we obtain the abstract differential equation 

d_ 

dt' 



: Q(t,s) = -£ (s-t)( A A y(t,s) + £ (s-t)(_ A 



-2W)£ () l«-l)[° Q I) foil -.s)Q(/..s) ->R(l.s)Q(l.s) (I.. 1.2) 



1 For a treatment without (Al) we construct £(t,s) = £o(t — s)Q(t,s) directly by the aid of formula 

(ESQ 
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with initial data Q(s,s) = I G C{H x H). Using that the operator norm satisfies 
s )||| < 2b(t) E L 1 [0,cxd) we see that Q(t,s) can be represented by the abstract 
Peano-Baker formula 

Q(t,s) = I + J2 K(h,s) K(t 2 ,s)... K(t k ,s)dt k ...dt 1 (1.13) 

k=l Js Js Js 

in terms of Bochner integrals. The L 1 -bound on the norm of 1Z implies a uniform bound 
on Q(i, s) and that 

lim Q(t,s) = Q(oo,s) (1.14) 

t^oo 

exists in C(H x i?) and satisfies 

/OO / /"OO \ 

6(r)dr exp f 2 y 6(0)d0j . (1.15) 

Furthermore, 

/OO / POO \ 

6(r)dr exp f 2 J b(9)d9j -» (1.16) 

as t — > oo implies invertibility of Q(oo,t) for sufficiently large while the propagation 
property Q(oo,s) = £q(s — t)Q(oo,t)£(t, s) extends this result to arbitrary s. 

Now the statement of Theorem 11.11 follows with JW+J* = lim^oo £q(— t)£(t, 0) = 
Q(oo,0) together with the fact that £o(t) is unitary. □ 

Remark 4. Because it was not necessary to use (Al), we can reverse the time direction in 
this statement and obtain the existence of a corresponding wave operator 

W_= lim J*£ (-t)£(t,0)J (1.17) 

t — > — OO 

when b G L 1 (M) and construct the scattering operator S = W+WI 1 . This operator maps 
incoming free waves (from t = — oo) to outgoing free waves (to t = oo), both parameterised 
by corresponding Cauchy data at time t = 0, which are connected by one solution of the 
dissipative equation. 

Remark 5. Theorem 1 1 . 1 1 remains true if we replace (jl.ip by 

u" + B(t)u + Au = 0, CCD) 
where B = B(t) G ^([0, oo), C(H)), cf. [6]. 
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2 Representation of solutions 



If b L 1 [0, oo) we will not apply (|1.13|) directly to represent the solutions of system (|1.8|) 
(although it is still valid). Our approach is based on a transformation of the system (|1.8p 
into a diagonal-dominated form. 

To make the calculations more transparent we use the spectral resolution of the operator 
A (or A, resp.) and reduce the analytic properties of the operator to algebraic properties 
of a corresponding function. Following [5l Theorem VIII. 4, p. 260] there exists a unitary 
operator S : H —> L 2 (X,dv), (X,dv) a suitable measure space, such that the operator 
A : H D T^{A) — > H is unitarily equivalent to 

(Au)(ti)=A(0u(0, u€V(A), (2.1) 

where -A(£) is a non-negative and ^-measurable function. Furthermore, 

V{A^)~L 2 {X,{l + A 2 \t))du), 7>0. (2.2) 

For the following the square root of A is of particular importance. We denote A(£) = 
■\/A(£). Then equation (II. ip is equivalent to 

d 2 t u{t , + 2b(t)d t u(t, + A 2 (£)u(t, = (2.3) 

and system (j 1 . 8 [) reads as 

(A(()u\ = I A({) \ (A(0u 



°> U» ) = {-HQ -W)) U» ) ■ {2A) 

Estimates in the energy space E correspond via H o J to estimates for this system in 
L 2 (X + ,dv) x L 2 (X,dv), where X + = {A(£) > 0}. 

We will solve this system following [8] and [9]. For a (later to be determined) number 
N > we decompose the extended phase space [0, oo) x X in two sets, the dissipative zone 

Z diss {N) := { (t, G [0, oo) x X \ A(£)(l + t) < N } (2.5) 

and the hyperbolic zone 

Z hyp (N) := { (t, <E [0, oo) x X \ A(^(l + t) > N }. (2.6) 

Let further % : X D {A(^) < A^} — > [0, oo) be implicitly defined by A(f )(1 + *{) = A^- 

In the first one we reformulate (|2.4p as system of Volterra integral equations and solve 
mainly by brute force, while in the second one we apply two steps of a diagonalisation 
procedure to extract the main terms of the representation of solutions and to get sharp 
estimates for the remainder terms. 
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2.1 Consideration in the dissipative zone 

We will prove estimates for the fundamental solution £(t,s,£) to (|2.4p . i.e. the matrix 
valued solution to that system with initial condition £(s,s,£) = I £ C 2x2 . Let for this 

A(t) = exp ( I b{T)dr\ (2.7) 



denote an auxiliary function related to the dissipation term b(t). We distinguish different 
cases related to the asymptotic behaviour of A(i). 

Case (CI): 71 := hm sup^^ tb(t) < 1/2. 

We denote by v(t,£) and w(t,£) the entries of one of the rows of £ (i, 0,£). Then these 
functions satisfy 

u(t,£) =7^ + A(0 / w(r,Odr, (2.8a) 
J o 

= a4) 772 " A(e) a4) / a2(tMt ' e)dr ' (2 ' 8b) 

where 77 = (771,772) = (1,0) or 77 = (0, 1) for the first and second column, respectively. 

Lemma 2.1. Assume (Al) and (A2). Then in the case (CI) the fundamental solution 
£ (t, 0,£) to (|2.4p satisfies for all (i,£) G Z(H SS (N) the point-wise estimat^ 



A*(t) VA- 7 (0 1, 
/or any 7 swc/t i/ta/j A 7 (^)A 2 (t) < 1. 

Example 2. If we consider the case 6(t) = /j/(l + i) with /j G [0, |) we obtain \(t) ~ i M and 
thus we can choose 7 = 2fi in the previous statement. 

Example 3. In general we can calculate the constant 7 as follows: Let 71 = limsup^^ to(i) 
and assume 71 < ^. Then the conditions of the above lemma are satisfied and X 2 (t) < i 2/i+e 
for any e > 0. Hence we choose 7 > 2/J. Especially, if tb(t) — > as i — ► 00 we have the 
above statement for any 7 > 0. Furthermore, if b G -L 1 [0, 00) the choice 7 = is possibleJl 

Proof, (of Lemma |2 . 1 j) We start by estimating the first column, i.e. 77 = (1,0). Plugging 
the second integral equation into the first one yields 

r* 1 



'* 1 



i-a 2 (0[ ^(o)v(e,o [ T^r-^dd, (2.10) 

Jo Je A 2 (t) 



For a matrix A we denote by \A\ the matrix consisting of the absolute values of the entries of A. 

_J I 1 

4(e+t) ~ (e + t)log(e + t) 



3 Otherwise, the choice 7 = 2/i is not always possible, as the example b(i) = 4 , e "V t , + , +t , 1[ | i / e+t ^ shows. 
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such that X 2 (t)A y (^)v(t, £) satisfies an Volterra integral equation with kernel k(t,9,£) = 
-A 2 (^)A 2 (t)] e *dr/A 2 (r) and source term h(t,g) = A 7 (£)A 2 (t) < 1. We represent its solu- 
tion by a Neumann series 



^ pt p 
r~; jo ■/ o 



00 /-t 

k(tk-i,t k , €)h(tk,€)dt k ■ ■ ■ dii. 



(2.11) 



From (CI) we know that t/X 2 (t) is monotone increasing for large t. Hence we conclude the 
kernel estimate 

pt pt ptf 

sup / sup \k(t,e,0\d9< sup A 2 (£)A 2 (t e ) / / Tir ^d9 
(t,OeZ dias (N) Jo o<t<^ (U)eZ dis3 {N) Jo Je X (rj 

= sup A 2 (OA 2 (t 5 ) f 6 -^—dr < A 2 (^f < 1 

(t,Z)eZ dUa (N) JO A t T ) 

which implies by standard arguments that | A 2 (i)A 7 (£)v (t, £)| < 1 uniformly on Z^ ss (iV). 
Then, the second integral equation implies 

\w(t,0\ < ^^A(0 j\x 2 (r)A^Ov(T,0\dr < ^^HOt < ( 2 - 12 ) 

by the definition of the zone. 

For the second column we use the same idea, plugging the second integral equation into 
the first one yields the new source term J* A(^)/A 2 (r)dr < 1/A 2 (i), such that \ 2 (t)v(t, 
is uniformly bounded on Zdi 8B (N). □ 

Case (C2): fj, := liminf^oo tb(t) > 1/2. 

In this case we obtain the decay rate 1 ft (which seems to be natural from the point of 
view of effective dissipation, cf. [10] and [TT]). 

Lemma 2.2. Assume (Al) and (A2). Then in the case (C2) the fundamental solution 
£(t,0,£) to (|2.4p satisfies for all (t,£) G Zdi S s{N) the point-wise estimate 

\£(t,0,0\ <T^-JVul IV (2-13) 



l + t \A- l (0 1, 

Proof. We proceed in a similar way like for Lemma 12.11 Plugging the second integral 
equation into the first one gives 

rt dr 



(1 + t)v(t, = (1 + t) m + (1 + t)A(0m [ 

Jo 



A 2 (r) 

■t 1 i + \2/ 



s 



Under condition (C2) we know that 1/A 2 (i) E L 1 (M_|_) and X 2 (t)/t is monotone increasing 
for large t. This implies again a kernel estimate for this Volterra integral equation 



rt rt _|_ £ _^2^^ 

SUp A 2 (£) / SUp / -i , a X 2(^\ dTd9 

-ZZM.JN) Jo o<i<tJe 1 + 



(t,S)ex d<ss (JV) ./o o<*<i e 

and we obtain for the first column that |A(£)(1 + t)v(t,£)\ < 1 is uniformly bounded on 
Zdiss{N) while for the second one |(1 + t)v(t,£)\ < 1 on Z^^N). The second integral 
equation yields corresponding bounds for w(t,£) from 

A «^i , IT7 dT ~ (1 + ()A(f) ~ 1 - <2 - 15) 

□ 

Remark 6. There remains a gap between the cases (CI) and (C2). In [7] the special case 
2b(t) = fi/(l + 1) for A = —A on L 2 (M. n ) was studied. The above exceptional value \i = 1 
is related to the occurrence of logarithmic terms in the representation of solutions. 

2.2 Consideration in the hyperbolic zone 

In Zhyp(N) we apply two steps of diagonalisation to system (12.4|> . In a first step we use 
to obtain 

In a second step we want to transform the second matrix without changing the structure 
of the first one. For this we set 

NM(t,0=i[ m (2-18) 

\2A(0 / 

B^(t,0 = d t N^{t,0 + b{t)N^. (2.19) 
Then we have by construction 



" iA<0 jAK) )^'(«)-^'(«)(- iAK) ^ -"*)(, (2.20) 
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such that with Ni(t, £) = I — N^(t, £) the operator identity 



holds. If the zone constant N is chosen sufficiently large, we have that det Ni(t,£) = 
1 — 6 2 (i)/(4A 2 (£)) > c > is uniformly bounded away from zero on Zy iyp {N). Then 
Nx (i 5 £) exists and N\(t,£) and A r 1 ~ 1 (t, ^) are both uniformly bounded on Zh yp (N). 
Setting R\(t,£) = — N^ 1 ^, ^-B*- 1 ) (t, £) we obtain the operator identity 



a - ( ' A(0 . A/t J+6(t)7 + i2i(t J (2.21) 



iA(0 

with remainder term i?i(t, £) subject to the point-wise estimate 

^ R ' (t -^^MtWTW (2 ' 22) 

following directly from Assumption (A2). Note, that we did not use Assumption (Al) for 
the treatment of the hyperbolic zone. 

We are now in a position to derive the main result of this section. 

Lemma 2.3. Assume (A2). Then the fundamental solution £(t,s,£) of (|2.4p can be rep- 
resented as 

£(t,s,0 = ^M~ 1 N^ 1 (t,0£o(t,s,OQi(t,s,ONi(s,OM (2.23) 
for t> s and (s,£) 6 Zh yp (N), where 



the function X(t) = exp y j Q 6(r)drJ describes the main influence of the dissipation 
b = b(t), 

the matrices Ni(t,£) and A^{" 1 (t,^) are uniformly bounded on Zh yp (N) tending on 
{A(£) > e} C X uniformly to the identity matrix I , 



the matrix £o(t,s,£) is given by 



-iA(0(ts) 



".-(/■■s.O ■ e <A( £ )(t- a ) ) . ( 2 - 24 ) 
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• and the matrix Qi(t,s,£) is uniformly bounded and invertible on Zh yp (N) tending on 
{A(£) > e} C X uniformly to the invertible matrix Qx(oo,s,£). 

Remark 7. The free propagator can be represented as £o(t—s, £) = M _1 <?o(i, s,£)M. Thus, 
Lemma [273] gives a relation between the free propagator £o(t — s, £) and £(t, s, £). 

Remark 8. If the operator A is boundedly invertible, the function A(£) is bounded from 
below and hence the matrices N\(t,£) and Qi(t,£) converge uniformly on X as t —* oo. 
Furthermore, the dissipative zone can be skipped in this case (because it is contained in 
a finite time strip). Thus Lemma 12.31 describes all essential properties of the fundamental 
solution in this case. 

Proof, (of Lemma 12. 3 j) The construction of the representation of solutions will be done in 
two steps. First, note that 

£o(t,s,Q = ^£ (t,s,0 (2.25) 
is the fundamental solution to the diagonal main part, 

dj (t, s, = f _lA ® . A(e) ) £ (t, s, + b(t)£o(t, s, 0, £o(s, s, = I. (2.26) 

Thus making the ansatz £o(t, s, £)Qi(i, s, £) for the fundamental solution to the transformed 
operator dt — diag (— «A(£), iA(£)) + b(t)I + R\(t, £), we obtain a system for Qi(t, s, £), 

d t Q 1 (t,s^)=TZ(t,s,OQi(t,s 1 0, Q 1 (s,8,$)=I (2.27) 

with coefficient matrix Tl(t,s,£) = £o(s, t, £)Ri(t, £)£o(t, s, £). Using that £o(t,s,£) is uni- 
tary, we see that TZ(t,s,£) satisfies the same estimates as £). 

This allows us to estimate in a second step the solution Qi(t, directly from the 
representation as Peano-Baker series 

Qi(t, s, = I + jr [ Kfas, f 1 K(t 2 , s, £) ■ • • r 1 n(t k , s, 0dt k ■ ■ ■ dti. (2.28) 

Thus, 

||Ci(t,fl,e)|| < exp (Y H^Cr.Olldr) < exp J* d-[_\ < 1 {2 29) 

uniformly in t > s and (s,£) G Zf iyp (N). Note that the inverse matrix satisfies a similar 
differential equation and therefore an analogue to series (12.281) . 
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It remains to show that Qi(t, s, £) converges uniformly on { A(£) >f}CXast-> oo. 
This follows by the Cauchy criterion applied to the series (|2.28|) or by the estimate 

/OO / ft 



< 1 f°° dr < 1 

~W)Jt (TTrj2 ~ e(TTt)' (2,30) 

The theorem is proven. □ 
2.3 Combination of results 

We combine the previously constructed representations of £(t, s, £) to get a representation 
for £ (t, 0, £) and the corresponding operator £(t, 0). We will obtain J*£(t, 0) J G and 



1 

Aft) 



|J^(t,0)J|||^<— (2.31) 



by the right choice of a smaller space E C E. Again J denotes the canonic identification 
E ~ (H GKerA) x H. 

If A(£) > N for all ( 6 I, we have to consider only the hyperbolic zone and £(t, 0,£) 
is given by Lemma 12.31 hence (I2.3ip holds true with E = E. We focus on the remaining 
case and consider £ with A(£) < N here. There it holds 

£(t,o,0 = £(t,k,0£(kAO (2.32) 

and we obtain corresponding estimates by Lemmata 12. II to 12.31 We distinguish some cases. 
Case (CI): Using Lemma [27T1 and [2T31 we conclude 



A(t) A^) VA" 7 (0 V ~ A(t) VA" 7 (0 1 

for A(£) < N,t>t{ and 7 with A 7 (£)A 2 (t ? ) < 1. For t < t$ the same estimate follows 
from Lemma 12.11 

To simplify notation we introduce [A(£)] = min(A(£), N). If we choose 7 = 1, we need 
for the data an estimate of the form ||[A(£)] _1 A(£)ui(£)||2 ~ ||(A(£))ui(£)|| 2 = ||«i||d(a)- 
Thus we obtain that the operator J*£(t,0)J restricted to T>(A) x H —* E satisfies the 
bound (J23H). 

Case (C2): Lemma 12.21 implies only the weaker decay rate 1/(1 + t). This can be com- 
pensated by assumptions on the data, which are valid near {A(£) = 0} C X. Using 

„ < 1 (2.34) 

1 + k (1 + ^)7+l ~ 1 > 
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for 7 > fj, — 1, fj, = limsup^^ tb(t), we conclude that 

for all A(£) < N and t > t%. For t < t^ Lemma 12.21 implies the same bound using 
A-r(OA(*)/(l + t)<A(t)/(l + *)T+ 1 <l. 

Thus if we define for 7 > the modified energy space EW to be the closure of (D(A) n 
i?(A 7 ) e Ker A) x (i?(A 7 ) Ker A) with respect to the norm 

\\( U1 ,u 2 )\\i {J) = WiHor^HOuiioWl^x,^ + wmr^MowUx^ v-m 

the operator J*£(t,0)J restricted to E^ — > E satisfies the bound (|2.3ip . 

Remark 9. Using [A(f)] _1 A(£) ~ (A(£)) we obtain especially £ (0) = (V(A) xH)Q(Ker A) 2 . 
So in the case (CI) we may use E^ for the data, while in (C2) we use E^ with 7 = 
max(7I + — 1,0), where 71 + denotes any number larger than Ji. 

Theorem 2.4. Let u be the solution to (jl.ip for data {u\, 112) G E^ forj = max(7i + — 1, 0). 
Then the estimate 

IIKnOlh^^HK,^)!!^) (2.37) 
is valid under assumptions (Al), (A 2) in the cases (CI) and (C2). 

Remark 10. We used conditions on the data related to {A(£) = 0}. In principle there arise 
two different cases. On the one hand, if i/{A(£) = 0} > it follows that G op{A) and 
is an eigenvalue of A. In this case we cut out the null-space of A and R{A) is dense in the 
ortho-complement of Ker A. Hence E^ is dense in E^ = (T> (A) Ker A) x (H Ker A), 
which is dense in the ortho-complement of {0} x Ker A in E. 

If v{A(^) = 0} = is a null-set, R(A) is dense in H and the result of Theorem 
holds on a dense subset of E^ = T>(A) x H. 



3 Scattering type theorems 

The main purpose of this note is to explain in which sense Theorem 12.41 is sharp. This 
sharpness is formulated as a generalisation of Theorem 11.11 with an additional energy decay 
function. 

Theorem 3.1. Under the assumptions of Theorem \2.J\ there exists a bounded operator 
W+ : £ (7) -> E, 7 = max(/I + - 1,0) such that for Cauchy data {u u u 2 ) G E^ of ([T7TT> 
and associated data (yi,V2) = W+ (1*1,1*2) to (|1.2|) the corresponding solutions u = u(t) and 
v = v(t) satisfy 

||A(t)(u,u')-(t;,t/)b-0 (3.1) 
as t — > 00. Furthermore, the operator W+ satisfies Ker W+ = {0}. 
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Proof. The proof is based on an explicit representation of the (modified) inverse wave 
operator W + . From Lemma Edl we know that the limit Qi(oo, if, £) = lim^oo Qi(i,tf,£) 
exists uniformly in {A(£) > e} for any e > 0. Hence, if we consider £o(— i, £)£(i, 0, £) on 
{A(£) > e} we obtain 

lim \(t)£o(-t,S)S(t,0,S) 

t—>oo 

= lim A(if)£ (-i, OM^N^it, £%{t, if, if , £)^i(if , 0^(^,0, 

= Hm A(t € )£b(-*, £)S>(* - if, O^Qifo % . 0^i(if , £)M£(if , 0, 
=A(if)£ (-if , 0M- l Qi{™, if , (if , OMSfc, 0, £) 
using A r -f 1 (i,^) — > I uniformly on {A(£) > e}. Denoting this limit by we see 

\w+(0\<HkMkA0\- (3-2) 

Now the estimates (|2T33]l and (|ll35j) on £(if,0,£) imply that W+{£) defines on (A(£) > 
0} = X \ {A(£) = 0} a multiplication operator with estimate 

IW+COI < !)[ A (0]" 7 (3-3) 

for 7 = max(/7 + — 1,0). This implies that the corresponding operator W+ = J*W+J 
resticted to E^' — > E is bounded. Furthermore, we obtain from the fact that the free 
propagator £o(t) is unitary for all t 

\\(ui,u 2 )\\ Eh) > \\\(t)£(t,0)(ui,u 2 ) - £ (t)W + ( Ul ,u 2 )\\ E 

= \\X(t)£ (-t)£{t,0)(u 1 ,u 2 )-W + {u 1 ,u 2 )\\ E -+0 (3.4) 

for all (ui,u 2 ) with suppf (ui, U2) C {A(£) > 0}. But this is a dense subset of E^ and by 
Banach-Steinhaus theorem this gives convergence for all (u\,u 2 ). 

On {A(£) > 0} the inverse wave operator W+ is represented as product of invertible 
matrices. Since {A(£) = 0} is excluded by the definition of E^\ the null-space Ker W + of 
W+ : E^ — > E is trivial and the theorem is proven. □ 

If we combine the previous theorem with the conservation of abstract energy for the 
free problem (|1.2|) . ||(v,v')||£; = II (vi,v 2 )\\e, we immediately obtain 

Corollary 3.2. For all data (ui,u 2 ) G E^> with 7 = max{7Z + — 1,0} the abstract energy 
of the solution to (jl.ip satisfies the two-sided estimate 

\\{u,u')\\ E ^^- (3.5) 

for all t > 0. 
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Note, that the definition of (formula (|2.36p ) implies _L KerA, such that 

corresponding eigenfunctions of A are excluded in the previous corollary. If the data 
belong to Ker A, the energy decays like 1/X 2 {t). 

If the operator A has a bounded inverse, i.e. A(£) > cq > uniformly on X, the 
dissipative zone Z,a ss (N) can be shrunken to Zdi ss (N)n{A(£) > cq}, such that 1+t < Ncq 1 
holds uniformly. In this case, the estimate of £(t,s,£) in the dissipative zone can be 
replaced by the trivial estimate \\£{t, s, £)|| < 1 following directly from system (|2,4p using 
the boundedness of \b(t)\ and A(£). Thus, in this case we can revert the time variable and 
consider the limits as t — > — oo under (A2) as well. Furthermore, E = E^ = T>(A) x H 
for all choices of 7 and Theorem 13.11 follows directly from Lemma 12.31 without relying on 
Banach-Steinhaus theorem. 

Theorem 3.3. Assume (Aw,w) > co\\w\\ for all w £ H. Assume further that the coeffi- 
cient b = b(t) satisfies (A2). Then there exist bounded and invertible operators W+ : E — > E 
and W- : E —> E, such that for all Cauchy data {u\,U2) G E and corresponding data 
(vf,vf) = W±(ui,U2) the estimate 

\\X(t)(u,u') - (v,v')\\ E < City'Wim,^)^ (3.6) 

holds true. The constant C depends on the coefficient b{t) and cq. 

Proof. The statement follows directly from Lemma 2.3 (for positive times and if we replace 
b{t) by —b{t) also for negative times) if we recall the following additional estimates (written 
for t > 0) 

\\N{\t,0-i\\ + ||Si(M,0 - Qi(oo,s,£)l| < A(e)(i + t) ~ TTt' (3 - 7) 

The first one follows from the properties of N^'(t,^), the last one restates (|2.30p . 

The operator W + is defined via the limit W+(^) as in the proof of Theorem 13.11 but 
the lower bound on A allows to skip the dissipative zone and to use te = instead. Then 
the above estimates imply for the difference 

\\\(t)£(t,o,0-£o(t)w + (0\\ 

= II M~ l N^\t, OM£ (t, M^Qiit, 0, $Ni (0, £)M 

S v ' S v ' 

=I+0(t~ 1 ) =(*)+0(t-i) 

- £ (t, M _1 Qi(oo, 0, ONi (0, OM \\ 

V v ' 

=(*) 

~ 1 + t 

and the desired statement follows. □ 
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Remark 11. In contrast to Theorem 13. II we can use Theorem 13.31 to obtain two-sided norm 
estimates for J*£(t, 0)J. It follows that 



J*£(t,0)J\\\ E ^E 



1 



(3.8) 



A(t) 



for all t. 

4 Applications 

We will sketch several applications to the previously given abstract results. The list is not 
complete, but it is intended to give an overview over the different cases and types of results. 

Example 4. As first example we consider A = —A on L 2 (M. n ). Then the abstract equation 
(|l.ip reads as wave equation with time-dependent dissipation. This case was treated in 
[8] and [9] by similar techniques. The abstract spaces reduce to Sobolev scale, we have 
P(A) = iI 1 (M") and E = H 1 ^ 1 ) x L 2 {W l ). 

In this special case Ker A = {0}. Thus, as special result, we obtain that for arbitrary 
data ui G ff x (R n ) and u 2 G L 2 (R n ) the hyperbolic energy of the corresponding solution 
u = u(t, x) satisfies 



as two-sided energy estimate. 

The same result is true for A = —A the Dirichlet- or Neumann-Laplacian on L 2 (Q) for 
an exterior domain £1 C R n with smooth boundary. 

Example 5. As second example we consider the Klein-Gordon type equation, where A = 
—A + 1 on L 2 (R n ). Then cr(A) = [1, oo) and A is boundedly invertible. Furthermore, 
X>(A) = H X (R) and E = iJ x (R n ) x L 2 (lR n ) and = E for all 7 > 0. 

In this case it is a simple consequence of the representation of Qi(t, s, £) that W + G C{E) 
is invertible in C(E). Especially (|4.ip is valid. 

The same result holds true in the case of A = —A the Dirichlet-Laplacian on L 2 (Q), Q 
an interior domain with smooth boundary. 

Example 6. Contrary to this for A = —A the Neumann-Laplacian on L 2 (£2), $7 bounded 
domain with smooth boundary, A has the non-trivial null-space Ker A = C consisting of 
constant functions. In this case the two sided estimate (|4.ip does not hold. Solving the 
corresponding ordinary differential equation we see that solutions to constant data behave 
like 1/X 2 {t) or b{t)/X 2 {t). 

Example 7. Let (aij(x))ij=i,...n G L°°(W n ,C nxn ) be a positive and self-adjoint matrix. 
Then we can consider the second order elliptic operator 



Vu(t,-)\\ 2 + \\u t (t,-)\\ 2 



1 



(4.1) 



X 2 (t) 




(4.2) 
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on L 2 (M n ) with domain T> {A) = H 2 (W l ). The corresponding abstract energy space has 
norm 



Example 8. The approach is not restricted to second order operators. We can also consider 
for example the damped plate equation with A = A 2 . In this case the abstract energy space 
is E = H 2 (W l ) x L 2 {R n ). 

Acknowledgements. The author thanks the referee for his many valuable comments to 
improve this paper. 
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